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1 Introduction 

In this paper we consider the Unary Subset-Sum problem which is defined as 
follows: Given integers mi, . . . ,m n and B (written in unary), we define the 
subset sum problem to be that of determining whether or not there exists an 
S C [n] so that X]ies m i — B (note that for this problem the mj are often 
assumed to be non-negative). Let C = \B\ + \ x i\ + This problem 

can be solved using a standard dynamic program using space 0(C) and time 
0(Cn). The dynamic program makes fundamental use of this large space and it 
is interesting to ask whether this requirement can be removed. Unary Subset- 
Sum has been studied in small-space models of computation as early as 1980 in 
[4], where they showed that it was in NL. Since then the problem was studied 
in [2] , where Cho and Huynh devised a complexity class between L and NL that 
contained Unary Subset-Sum as supporting evidence that it is not iVX-complete. 
This problem was listed again in pQ claiming it to be an open problem as to 
whether or not it is in L. It was recently shown in [3] that this problem was in 
Logspace. Unfortunately, their algorithm is somewhat complicated. We provide 
a simple algorithm solving this problem in Logspace. 

2 Our Algorithm 

The basic idea of our algorithm will be to make use of the generating function 
n™=i(l + x" li ) — J2sc[n] x^its m i to compute the number of solutions to our 
problem modulo p for a number of different primes p (we show how to do this 
in Lemma [l}. Pseudocode for our algorithm is follows: 
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c := 

p := NcxtPrime(C) 
While(c < n) 

If El'Ji *- B nr=i(l + * mi ) # (mod p) 
Return True 

c := c+ Llog 2 (p)J 

p := NcxtPrime(p) 
End While 
Return False 



2.1 Space Complexity 

There are several things that must be noted to show that this algorithm runs 
in logspace. First we claim that p is never more than polynomial in size. This 
is because standard facts about prime numbers imply that there are at least n 
primes between C and poly(C, n), and each of these primes causes c to increase 
by at least 1. We also note that Y^x=i X ~ B Y\a=i(^ + x mi ) can be computed 
modulo p in Logspace. This is because we can just keep track of the value 
of x and the current running total (modulo p) along with the space necessary 
to compute the next term. The product is computed again by keeping track 
of i and the current running product (modulo p) and whatever is necessary 
to compute the next term. The exponents are computed in the obvious way. 
Finally primality testing of poly-sized numbers can be done by repeated trial 
divisions in Logspace, and hence the NextPrime function can also be computed 
in Logspace. 

2.2 Correctness 

We now have to prove correctness of the algorithm. Let A be the number of 
subsets S C [n] so that Eies TO i = B. 

Lemma 1. For p a prime number, p > C. Then 

p-l n 

J2 X ~ B H( 1 + x " li ) = ~ A (modp). 

x—l i—1 

Where again A is the number of subsets S C [n] so that J2ieS m i = B . 
Proof. Note that 

n 

X- B \\{l+X m >)= xZiesm-B. 
t=l SC[n] 

The idea of our proof will be to interchange the order of summation and show 
that the terms for which Eies m « B cancel out. 



2 



Notice that each exponent in this sum has absolute value less than p — 1. 
Interchanging the sums on the right hand side, we find that 



p — 1 n p— 1 



]>>-*n(i+- mi )= e E^ esi 

x=l i=l SC[n] i=l 



We note that: 



£V (modp),^ 1 ( mod P - !) 



x=l 



If fc is a multiple of p — 1, then all terms in the sum are 1 modulo p and the 
result follows. Otherwise, we let g be a primitive root mod p and note that 
instead of summing over x — 1 to p — 1 we may sum over x = g £ for t = to 
p - 2. Then 

Hence 

p— 1 n p— 1 

£V B I](i E E 2 ^" 1 " 5 ^ E - L 

x=l i=l Sc[n]i=l SC[n] 

EieS x i=B (mod p— 1) 

Since p— 1 is larger than C, X^ies x i = B (mod p— 1) if and only if XieS = 
Hence this sum contributes -1 for each such S and so the final sum is —A. □ 

We are now ready to prove correctness. If Ylx=i X ~ B nr=i(l + x mi ) ^ 
(mod p) for some p > C, then by our Lemma, this means that A ^ (mod p) . 
In particular, this means that A ^ 0, and that therefore there is some such 
S. Consider an integer d which is equal to the product of the primes p that 
have been checked so far. Then d is a product of distinct primes p so that 
—A = Xx=i X ~ B n"=i(l + x mi ) = (mod p). Therefore d\A. Furthermore it 
is the case that d > 2 C . It is clear from the definition of A that < A < 2™. 
Therefore if c > n, d > 2™ and d\A, which implies that A = 0, and that therefore 
there are no solutions. Hence our algorithm always outputs correctly. 



3 Extensions 

There are some relatively simple extensions of this algorithm. For one thing, 
our algorithm does more than tell us whether or not A is equal to 0, but also 
tells us congruential information about A. We can in fact obtain more refined 
congrucntial information than is apparent from our Lemma. We can also use 
this along with the Chinese Remainder Theorem to compute a numerical ap- 
proximation of A. Finally a slight generalization of these techniques allows us 
to work with to, vector valued rather than integer-valued. 



3 



3.1 Computing Congruences 



We show above how to compute A modulo p for p a prime larger than C. But 
in fact if p is any prime and k > 1 any integer, A can be computed modulo p k 
in 0(\og(p k )) space. 

If p > C, then we have that 

1 p-i n 

1 X=l 1=1 

On the other hand if p < C, the above expression will only count the number 
of subsets that give the correct sum modulo p — 1. We can fix this by letting 
q = p l for some integer I so that q > C. Then for the same reasons that the 
above is true, it will be the case that 

n 

A= — - x- B Y[(l + x" H ) (modp). 

q xe¥* i=l 

Where ¥ q is the finite field of order q. 

If we have k > 1 and p > C we note that again for the same reasons 

^"^E x- B f[(l + x^) (mod/). 

xEfJ.p-1 i=l 

Where [i p -\ is the set of (p — l) st roots of unity in 7Ljp k . This computation can 
be performed without difficulty in "L/p k . We again run into difficulty if p < C. 
We can let r = (p — \)p l for some integer I so that r > C. It will then be the 
case that 

n 

rA=J2 x- B + x mi ) (mod p k+e ). 

i£(i r i=l 

The right hand side can easily be computed in and dividing by r gives 

A (mod p k ). 

3.2 Approximating the Number of Solutions 

It is also possible in Logspace to approximate the number of solutions, A, com- 
puting logarithmically many significant bits. This can be done using the Chinese 
Remainder Theorem. Suppose that pi, . . . ,pt are distinct primes. By the above 
we can compute A modulo pi for each i. Let N — Yl^iPii an d Ni = The 
Chinese Remainder Theorem tells us that 

k 

A = ^2,N i {A (mod Pl )) (Nr 1 (mod ft)) (mod N). 
»=i 
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Or in other words, 



N=J2[-)( A (mod (TVri ( mo d Pi )) (modi). 

Now we can compute A modulo pi by the above. We can also compute A^ 1 = 
Yij^iPj 1 (mod pi). Hence we can compute each term in the sum to logarith- 
mically many bits. Hence in logspace we can compute 

| ( modl ) 

to logarithmically many bits of precision. If 2A > N > A, this allows us to 
compute logarithmically many significant bits of A. We can find such an N by 
starting with an N > 2™ > A and repeatedly trying N at least half as big as 
the previous N until N < 2A (we can find our next N by either removing the 
prime 2 from N or replacing the smallest prime dividing N by one at least half 
as big (which exists by Bertrand's postulate)). 

It should also be noted that this ability to approximately count solutions 
in Logspace allows us to approximately uniformly sample from the space of 
solutions in Randomized Logspace. This is done by deciding whether or not 
each element is in S one-by-one and putting it in with probability nearly equal 
to the proportion of the remaining solutions that have that element in S. 



3.3 Vector- Valued Inputs 

We consider the slightly modified subset sum problem where now rrii and B 
lie in Z fe , and again we wish to determine whether or not there exists and S 
so that X^ieS mi = ^- If we let C be one more than the sum of the absolute 
values of the coefficients of the rm plus the absolute values of the coefficients 
of B, a slight modification of our algorithm allows us to solve this problem 
in 0(fclog(C)) space and time (in particular if k = 0(1), this runs in 

0(log(C)) space and C° {1) time). 

There are two ways to do this. One is simply to treat our vectors as base 
C-expansions of integers and reduce this to our previous algorithm. Another 
technique involves a slight generalization of our Lemma. In either case we let 
mi = (mi,i, . . .,mi,fe), B = (B 1 , . . .,B k ). 

For the first algorithm, we let = X^j=i ^~ lm %,j an d B' — X^j=i 
We claim that for any S C [n] that ^2 ieS mi — B if and only if J2ies m 'i = 
thus reducing this to an instance of our original problem. The claim holds 
because 

£ mj - B' = £ (j2 mtj Bj) = £ ^ V 

ies 3=1 Vies / j=i 

Since the e 3 are all integers of absolute value less than C, this sum is if and 
only if, each of the ej are 0. Hence J2i<=s m i = B if and only if J2ies m i = B' ■ 
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Another way to do this is by generalizing our Lemma. In particular it can 
be shown using similar techniques that if A is the number of subsets S that 
work, and if p is a prime bigger than C that 



-A = 



xi,...,Xk — l \i— 1 



p-1 / k 



e n 



-Bi 



) 



( 



i=l \ j=l 



n I k 

n i+ ik 




(mod p). 



Given this, there is a natural generalization of our algorithm. 

It should also be noted that both of these techniques allow us to use the 
above-stated generalizations to our algorithm in the vector-valued context. 

This generalization also allows us to solve some related problems, such as 
the Unary 0-1 Knapsack problem. This problem is defined as follows: You are 
given a list of integer weights w\, . . . , w n , a list of integer values, v\, . . . , v n , and 
an integer bound B. The objective is to find a subset S C [n] so that^ ieS m is 
as large as possible subject to the restriction that X^es Wi — We do this by 
determining all possible pairs of (J2ies ~52ieS Vi ) ^ a PPly m g our algorithm 
to nii — {wi,Vi) and B — (w,v) for all \w\ < X)"=i 1^1' \ v \ — 1^1- ^f the 

pairs (w, v) for which there is a solution, we keep track of the largest v that 
corresponds to a w < B. From this pair it is also not hard to use our algorithm 
to find a subset S which achieves this bound. 
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